Blood RNA alternative splicing events as diagnostic biomarkers for infectious disease

Summary Assays detecting blood transcriptome changes are studied for infectious disease diagnosis. Blood-based RNA alternative splicing (AS) events, which have not been well characterized in pathogen infection, have potential normalization and assay platform stability advantages over gene expression for diagnosis. Here, we present a computational framework for developing AS diagnostic biomarkers. Leveraging a large prospective cohort of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and whole-blood RNA sequencing (RNA-seq) data, we identify a major functional AS program switch upon viral infection. Using an independent cohort, we demonstrate the improved accuracy of AS biomarkers for SARS-CoV-2 diagnosis compared with six reported transcriptome signatures. We then optimize a subset of AS-based biomarkers to develop microfluidic PCR diagnostic assays. This assay achieves nearly perfect test accuracy (61/62 = 98.4%) using a naive principal component classifier, significantly more accurate than a gene expression PCR assay in the same cohort. Therefore, our RNA splicing computational framework enables a promising avenue for host-response diagnosis of infection. Graphical abstract Highlights • We present a computational framework for alternative splicing (AS) diagnostic markers• Our AS biomarkers outperform gene-expression biomarkers in COVID-19 detection• Microfluidic PCR diagnostic assay of AS biomarkers achieves greater than 98% accuracy• We interpret the biological importance of identified AS biomarkers Motivation Host-based response assays (HRAs) can often diagnose infectious disease earlier and more precisely than pathogen-based tests. However, the role of RNA alternative splicing (AS) in HRAs remains unexplored, as existing HRAs are restricted to gene expression signatures. We report a computational framework for the identification, optimization, and evaluation of blood AS-based diagnostic assay development for infectious disease. Using SARS-CoV-2 infection as a case study, we demonstrate the improved accuracy of AS biomarkers for COVID-19 diagnosis when compared against six reported transcriptome signatures and when implemented as a microfluidic PCR diagnostic assay. Host-based response assays can diagnose infectious disease earlier and more precisely than pathogen-based tests. However, the role of RNA alternative splicing (AS) remains unexplored. Zhang et al. present a computational framework for AS diagnostic biomarkers. Using SARS-CoV-2 as a case study, they demonstrate the improved accuracy of AS biomarkers for COVID-19 diagnosis.

Blood RNA alternative splicing events as diagnostic biomarkers for infectious disease.

Assays detecting blood transcriptome changes are studied for infectious disease diagnosis. Blood-based RNA alternative splicing (AS) events, which have not been well characterized in pathogen infection, have potential normalization and assay platform stability advantages over gene expression for diagnosis. Here, we present a computational framework for developing AS diagnostic biomarkers. Leveraging a large prospective cohort of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and whole-blood RNA sequencing (RNA-seq) data, we identify a major functional AS program switch upon viral infection. Using an independent cohort, we demonstrate the improved accuracy of AS biomarkers for SARS-CoV-2 diagnosis compared with six reported transcriptome signatures. We then optimize a subset of AS-based biomarkers to develop microfluidic PCR diagnostic assays. This assay achieves nearly perfect test accuracy (61/62 = 98.4%) using a naive principal component classifier, significantly more accurate than a gene expression PCR assay in the same cohort. Therefore, our RNA splicing computational framework enables a promising avenue for host-response diagnosis of infection.

Multi-objective optimization identifies a specific and interpretable COVID-19 host response signature.

The identification of a COVID-19 host response signature in blood can increase the understanding of SARS-CoV-2 pathogenesis and improve diagnostic tools. Applying a multi-objective optimization framework to both massive public and new multi-omics data, we identified a COVID-19 signature regulated at both transcriptional and epigenetic levels. We validated the signature's robustness in multiple independent COVID-19 cohorts. Using public data from 8,630 subjects and 53 conditions, we demonstrated no cross-reactivity with other viral and bacterial infections, COVID-19 comorbidities, or confounders. In contrast, previously reported COVID-19 signatures were associated with significant cross-reactivity. The signature's interpretation, based on cell-type deconvolution and single-cell data analysis, revealed prominent yet complementary roles for plasmablasts and memory T cells. Although the signal from plasmablasts mediated COVID-19 detection, the signal from memory T cells controlled against cross-reactivity with other viral infections. This framework identified a robust, interpretable COVID-19 signature and is broadly applicable in other disease contexts. A record of this paper's transparent peer review process is included in the supplemental information.

Systematic comparison of published host gene expression signatures for bacterial/viral discrimination.

BACKGROUND: Measuring host gene expression is a promising diagnostic strategy to discriminate bacterial and viral infections. Multiple signatures of varying size, complexity, and target populations have been described. However, there is little information to indicate how the performance of various published signatures compare to one another. METHODS: This systematic comparison of host gene expression signatures evaluated the performance of 28 signatures, validating them in 4589 subjects from 51 publicly available datasets. Thirteen COVID-specific datasets with 1416 subjects were included in a separate analysis. Individual signature performance was evaluated using the area under the receiving operating characteristic curve (AUC) value. Overall signature performance was evaluated using median AUCs and accuracies. RESULTS: Signature performance varied widely, with median AUCs ranging from 0.55 to 0.96 for bacterial classification and 0.69-0.97 for viral classification. Signature size varied (1-398 genes), with smaller signatures generally performing more poorly (P < 0.04). Viral infection was easier to diagnose than bacterial infection (84% vs. 79% overall accuracy, respectively; P < .001). Host gene expression classifiers performed more poorly in some pediatric populations (3 months-1 year and 2-11 years) compared to the adult population for both bacterial infection (73% and 70% vs. 82%, respectively; P < .001) and viral infection (80% and 79% vs. 88%, respectively; P < .001). We did not observe classification differences based on illness severity as defined by ICU admission for bacterial or viral infections. The median AUC across all signatures for COVID-19 classification was 0.80 compared to 0.83 for viral classification in the same datasets. CONCLUSIONS: In this systematic comparison of 28 host gene expression signatures, we observed differences based on a signature's size and characteristics of the validation population, including age and infection type. However, populations used for signature discovery did not impact performance, underscoring the redundancy among many of these signatures. Furthermore, differential performance in specific populations may only be observable through this type of large-scale validation.

SARS-CoV-2 reinfection across a spectrum of immunological states.

Purpose: Several cases of symptomatic reinfection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) after full recovery from a prior episode have been reported. As reinfection has become an increasingly common phenomenon, an improved understanding of the risk factors for reinfection and the character and duration of the serological responses to infection and vaccination is critical for managing the coronavirus disease 2019 (COVID-19) pandemic. Methods: We described four cases of SARS-CoV-2 reinfection in individuals representing a spectrum of healthy and immunocompromised states, including (1) a healthy 41-year-old pediatrician, (2) an immunocompromised 31-year-old with granulomatosis with polyangiitis, (3) a healthy 26-year-old pregnant woman, and (4) a 50-year-old with hypertension and hyperlipidemia. We performed confirmatory quantitative reverse transcription-polymerase chain reaction and qualitative immunoglobulin M and quantitative IgG testing on all available patient samples to confirm the presence of infection and serological response to infection. Results: Our analysis showed that patients 1 and 2, a healthy and an immunocompromised patient, both failed to mount a robust serologic response to the initial infection. In contrast, patients 3 and 4, with minimal comorbid disease, both mounted a strong serological response to their initial infection, but were still susceptible to reinfection. Conclusion: Repeat episodes of COVID-19 are capable of occurring in patients regardless of the presence of known risk factors for infection or level of serological response to infection, although this did not trigger critical illness in any instance.

Differential chromatin accessibility in peripheral blood mononuclear cells underlies COVID-19 disease severity prior to seroconversion.

SARS-CoV-2 infection triggers profound and variable immune responses in human hosts. Chromatin remodeling has been observed in individuals severely ill or convalescing with COVID-19, but chromatin remodeling early in disease prior to anti-spike protein IgG seroconversion has not been defined. We performed the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) and RNA-seq on peripheral blood mononuclear cells (PBMCs) from outpatients with mild or moderate symptom severity at different stages of clinical illness. Early in the disease course prior to IgG seroconversion, modifications in chromatin accessibility associated with mild or moderate symptoms were already robust and included severity-associated changes in accessibility of genes in interleukin signaling, regulation of cell differentiation and cell morphology. Furthermore, single-cell analyses revealed evolution of the chromatin accessibility landscape and transcription factor motif accessibility for individual PBMC cell types over time. The most extensive remodeling occurred in CD14+ monocytes, where sub-populations with distinct chromatin accessibility profiles were observed prior to seroconversion. Mild symptom severity was marked by upregulation of classical antiviral pathways, including those regulating IRF1 and IRF7, whereas in moderate disease, these classical antiviral signals diminished, suggesting dysregulated and less effective responses. Together, these observations offer novel insight into the epigenome of early mild SARS-CoV-2 infection and suggest that detection of chromatin remodeling in early disease may offer promise for a new class of diagnostic tools for COVID-19.

Prospective Validation of a Rapid Host Gene Expression Test to Discriminate Bacterial From Viral Respiratory Infection.

Importance: Bacterial and viral causes of acute respiratory illness (ARI) are difficult to clinically distinguish, resulting in the inappropriate use of antibacterial therapy. The use of a host gene expression-based test that is able to discriminate bacterial from viral infection in less than 1 hour may improve care and antimicrobial stewardship. Objective: To validate the host response bacterial/viral (HR-B/V) test and assess its ability to accurately differentiate bacterial from viral infection among patients with ARI. Design, Setting, and Participants: This prospective multicenter diagnostic study enrolled 755 children and adults with febrile ARI of 7 or fewer days' duration from 10 US emergency departments. Participants were enrolled from October 3, 2014, to September 1, 2019, followed by additional enrollment of patients with COVID-19 from March 20 to December 3, 2020. Clinical adjudication of enrolled participants identified 616 individuals as having bacterial or viral infection. The primary analysis cohort included 334 participants with high-confidence reference adjudications (based on adjudicator concordance and the presence of an identified pathogen confirmed by microbiological testing). A secondary analysis of the entire cohort of 616 participants included cases with low-confidence reference adjudications (based on adjudicator discordance or the absence of an identified pathogen in microbiological testing). Thirty-three participants with COVID-19 were included post hoc. Interventions: The HR-B/V test quantified the expression of 45 host messenger RNAs in approximately 45 minutes to derive a probability of bacterial infection. Main Outcomes and Measures: Performance characteristics for the HR-B/V test compared with clinical adjudication were reported as either bacterial or viral infection or categorized into 4 likelihood groups (viral very likely [probability score <0.19], viral likely [probability score of 0.19-0.40], bacterial likely [probability score of 0.41-0.73], and bacterial very likely [probability score >0.73]) and compared with procalcitonin measurement. Results: Among 755 enrolled participants, the median age was 26 years (IQR, 16-52 years); 360 participants (47.7%) were female, and 395 (52.3%) were male. A total of 13 participants (1.7%) were American Indian, 13 (1.7%) were Asian, 368 (48.7%) were Black, 131 (17.4%) were Hispanic, 3 (0.4%) were Native Hawaiian or Pacific Islander, 297 (39.3%) were White, and 60 (7.9%) were of unspecified race and/or ethnicity. In the primary analysis involving 334 participants, the HR-B/V test had sensitivity of 89.8% (95% CI, 77.8%-96.2%), specificity of 82.1% (95% CI, 77.4%-86.6%), and a negative predictive value (NPV) of 97.9% (95% CI, 95.3%-99.1%) for bacterial infection. In comparison, the sensitivity of procalcitonin measurement was 28.6% (95% CI, 16.2%-40.9%; P < .001), the specificity was 87.0% (95% CI, 82.7%-90.7%; P = .006), and the NPV was 87.6% (95% CI, 85.5%-89.5%; P < .001). When stratified into likelihood groups, the HR-B/V test had an NPV of 98.9% (95% CI, 96.1%-100%) for bacterial infection in the viral very likely group and a positive predictive value of 63.4% (95% CI, 47.2%-77.9%) for bacterial infection in the bacterial very likely group. The HR-B/V test correctly identified 30 of 33 participants (90.9%) with acute COVID-19 as having a viral infection. Conclusions and Relevance: In this study, the HR-B/V test accurately discriminated bacterial from viral infection among patients with febrile ARI and was superior to procalcitonin measurement. The findings suggest that an accurate point-of-need host response test with high NPV may offer an opportunity to improve antibiotic stewardship and patient outcomes.

Assessment of the Feasibility of Using Noninvasive Wearable Biometric Monitoring Sensors to Detect Influenza and the Common Cold Before Symptom Onset.

Importance: Currently, there are no presymptomatic screening methods to identify individuals infected with a respiratory virus to prevent disease spread and to predict their trajectory for resource allocation. Objective: To evaluate the feasibility of using noninvasive, wrist-worn wearable biometric monitoring sensors to detect presymptomatic viral infection after exposure and predict infection severity in patients exposed to H1N1 influenza or human rhinovirus. Design, Setting, and Participants: The cohort H1N1 viral challenge study was conducted during 2018; data were collected from September 11, 2017, to May 4, 2018. The cohort rhinovirus challenge study was conducted during 2015; data were collected from September 14 to 21, 2015. A total of 39 adult participants were recruited for the H1N1 challenge study, and 24 adult participants were recruited for the rhinovirus challenge study. Exclusion criteria for both challenges included chronic respiratory illness and high levels of serum antibodies. Participants in the H1N1 challenge study were isolated in a clinic for a minimum of 8 days after inoculation. The rhinovirus challenge took place on a college campus, and participants were not isolated. Exposures: Participants in the H1N1 challenge study were inoculated via intranasal drops of diluted influenza A/California/03/09 (H1N1) virus with a mean count of 106 using the median tissue culture infectious dose (TCID50) assay. Participants in the rhinovirus challenge study were inoculated via intranasal drops of diluted human rhinovirus strain type 16 with a count of 100 using the TCID50 assay. Main Outcomes and Measures: The primary outcome measures included cross-validated performance metrics of random forest models to screen for presymptomatic infection and predict infection severity, including accuracy, precision, sensitivity, specificity, F1 score, and area under the receiver operating characteristic curve (AUC). Results: A total of 31 participants with H1N1 (24 men [77.4%]; mean [SD] age, 34.7 [12.3] years) and 18 participants with rhinovirus (11 men [61.1%]; mean [SD] age, 21.7 [3.1] years) were included in the analysis after data preprocessing. Separate H1N1 and rhinovirus detection models, using only data on wearble devices as input, were able to distinguish between infection and noninfection with accuracies of up to 92% for H1N1 (90% precision, 90% sensitivity, 93% specificity, and 90% F1 score, 0.85 [95% CI, 0.70-1.00] AUC) and 88% for rhinovirus (100% precision, 78% sensitivity, 100% specificity, 88% F1 score, and 0.96 [95% CI, 0.85-1.00] AUC). The infection severity prediction model was able to distinguish between mild and moderate infection 24 hours prior to symptom onset with an accuracy of 90% for H1N1 (88% precision, 88% sensitivity, 92% specificity, 88% F1 score, and 0.88 [95% CI, 0.72-1.00] AUC) and 89% for rhinovirus (100% precision, 75% sensitivity, 100% specificity, 86% F1 score, and 0.95 [95% CI, 0.79-1.00] AUC). Conclusions and Relevance: This cohort study suggests that the use of a noninvasive, wrist-worn wearable device to predict an individual's response to viral exposure prior to symptoms is feasible. Harnessing this technology would support early interventions to limit presymptomatic spread of viral respiratory infections, which is timely in the era of COVID-19.

Mucosal-associated invariant T cell responses differ by sex in COVID-19.

BACKGROUND: Sexual dimorphisms in immune responses contribute to coronavirus disease 2019 (COVID-19) outcomes, but the mechanisms governing this disparity remain incompletely understood. METHODS: We carried out sex-balanced sampling of peripheral blood mononuclear cells from hospitalized and non-hospitalized individuals with confirmed COVID-19, uninfected close contacts, and healthy control individuals for 36-color flow cytometry and single-cell RNA sequencing. FINDINGS: Our results revealed a pronounced reduction of circulating mucosal-associated invariant T (MAIT) cells in infected females. Integration of published COVID-19 airway tissue datasets suggests that this reduction represented a major wave of MAIT cell extravasation during early infection in females. Moreover, MAIT cells from females possessed an immunologically active gene signature, whereas cells from males were pro-apoptotic. CONCLUSIONS: Our findings uncover a female-specific protective MAIT cell profile, potentially shedding light on reduced COVID-19 susceptibility in females. FUNDING: This work was supported by NIH/NIAID (U01AI066569 and UM1AI104681), the Defense Advanced Projects Agency (DARPA; N66001-09-C-2082 and HR0011-17-2-0069), the Veterans Affairs Health System, and Virology Quality Assurance (VQA; 75N93019C00015). The content is solely the responsibility of the authors and does not necessarily represent the official view of the National Institutes of Health. COVID-19 samples were processed under Biosafety level 2 (BSL-2) with aerosol management enhancement or BSL-3 in the Duke Regional Biocontainment Laboratory, which received partial support for construction from NIH/NIAID (UC6AI058607).

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility.

BACKGROUND: While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility. RESULTS: Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb)) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs associated with both severe COVID-19 and other human traits demonstrated colocalization of the GWAS signal at the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN). This finding points to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity. CONCLUSIONS: Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches. The iCPAGdb web portal is accessible at http://cpag.oit.duke.edu and the software code at https://github.com/tbalmat/iCPAGdb .

Dysregulated transcriptional responses to SARS-CoV-2 in the periphery.

SARS-CoV-2 infection has been shown to trigger a wide spectrum of immune responses and clinical manifestations in human hosts. Here, we sought to elucidate novel aspects of the host response to SARS-CoV-2 infection through RNA sequencing of peripheral blood samples from 46 subjects with COVID-19 and directly comparing them to subjects with seasonal coronavirus, influenza, bacterial pneumonia, and healthy controls. Early SARS-CoV-2 infection triggers a powerful transcriptomic response in peripheral blood with conserved components that are heavily interferon-driven but also marked by indicators of early B-cell activation and antibody production. Interferon responses during SARS-CoV-2 infection demonstrate unique patterns of dysregulated expression compared to other infectious and healthy states. Heterogeneous activation of coagulation and fibrinolytic pathways are present in early COVID-19, as are IL1 and JAK/STAT signaling pathways, which persist into late disease. Classifiers based on differentially expressed genes accurately distinguished SARS-CoV-2 infection from other acute illnesses (auROC 0.95 [95% CI 0.92-0.98]). The transcriptome in peripheral blood reveals both diverse and conserved components of the immune response in COVID-19 and provides for potential biomarker-based approaches to diagnosis.